Is the Sgr dSph a dark matter dominated system? 
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Abstract 

We study the evolution of possible progenitors of Sgr dSph using several numerical N-body simulations of different 
dwarf spheroidal galaxies both with and without dark matter, as they orbit the Milky Way. The barionic and dark 
components of the dwarfs were made obeying a Plummer and NFW potentials of 10 6 particles respectively. The 
Milky Way was modeled like a tree-component rigid potential and the simulations were performed using a modified 
Gadget-2 code. We found that none of the simulated galaxies without dark matter reproduced the physical properties 
observed in Sgr dSph, suggesting that, at the beginning of its evolution, Sgr dSph might have been immersed in a dark 
matter halo. 

The simulations of progenitors immersed in dark matter halos suggest that Sgr dSph at its beginning might have been 
an extended system, i.e. its Plummer radius could have had a value approximated to 1.2 kpc or higher; furthermore, 
this galaxy could have been immersed in a dark halo with a mass higher than 10 8 M Q . These results are important for 
the construction of a model of the formation of Sgr dSph. 
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1. Introduction 

Since its discovery in 1994 by Ibata, Gilmore and Irwin [ 19], the Sgr dSph dwarf galaxy has been object of febrile 
research; but only recently, thanks to the new data provided by 2MASS and SDSS, its tidal tails have been studied in 
great detail Il43l IT61 [32l [36l n~3l l4l [5X1 1501 . The Sgr dSph is very important because it gives us a better understanding 
about the formation process of the Milky Way and many clues about the shape of its dark matter halo. 

Sgr dSph tidal tails are the biggest contributions to the galactic halo lfT4l at a radius above 50 kpc 11291 . Many 
studies have found tidally stripped material associated to this dwarf with an age between 3 and 6.5 Gyrs [4|; others, 
suggest that the age of the older debris is between 2 and 4 Gyrs [ 17 1; furthermore, arc substructures have been found 
through the galactic south hemisphere at a distance of approximately 40 kpc ll36l : as well as filaments at a distance of 
45 kpc E9l|38l; 46 ± 12 kpc ED and 62 ± 6 kpc [38 1 to the galactic center. These data suggest that Sgr dSph has been 
disrupted since about 6 Gyrs. 

The nearest Sgr dSph tidal tails have been found at 18 kpc from the galactic center. These debris belong to fil- 
aments that emerged during the last pericenter passage of the satellite, approximately 0.75 Gyrs ago 1 10|. The most 
recent discoveries about the structure of the tidal tails of the Sgr dSph, were done by Belokurov et al. [4] and Yanny 
et al. [51 1 who found two branches of the Sgr dSph stream -the leading and trailing tidal tails- in an area around the 
north galactic cap. 
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Table 1: Physical properties observed in Sgr dSph by several 
authors. The nomenclature is referred to: * I41I , ** 1211 . *** 

UH ED 

H/2 0.55 kpc*** 

o- 11.4±7km/s** 

yUo.v 25.4 mag/arcsec 2 *** 
2.47 x 10 6 L Q /kpc 2 

M/L 10* 
14-16 + 



Sgr dSph main body is located at 16 ± 2 kpc from the galactic center ll2Tll24ll20ll22l[T5ll . at a distance of 24 ± 2 
kpc from the sun |27][3]|2l] with / = 5°, b = -14.5° EQIU. This dwarf covers a vast region on the sky, 15° x 7° lfl4l . 
Numerical simulations suggest that it is not symmetric; that is, the main body of Sgr dSph is a prolate object with axis 
ratios * 3:1:1 ||20l|22l. 

The central velocity dispersion of the main body of Sgr dSph is 11.4 + 0.7 km/s l20l |2T1 . Kinematical studies 
suggest that Sgr dSph moves to the north with transversal velocity of 250 ± 90 km/s Il25ll2l1l . The radial velocity of 
the satellite is 171 km/s l25l l28l fT5l 1571 l32l l33l . 

Studies realized by Law et al. [32 1 and Majewski et al. [37] indicate that the mass bound to the main body of Sgr 
dSph is a: 3 x 10 8 M Q . A better estimate was done by Law and Majewski |[35], who conclude that M = 2.5+j 3 x 10 8 
M Q ; with this range of masses, the absolute visual magnitude of the satellite is M v = -13.64; its total luminosity, 
Lj = 2.4 x 1O 7 L and its mass to light ratio, M/L = 10^ M G /L Q . Table [l] shows some physical properties of this 
dwarf observed by other authors. r;/2 is the half light radius; cr is the central line of sight velocity dispersion; jUq is 
the central surface brightness and M/L is the mass to light ratio. 

The aim of this paper is to show possible features of the progenitor of Sgr dSph and an estimate of the mass 
of its initial dark matter halo using N-body simulations. We present our model in section [2] the results obtained 
by simulating several possible progenitors of Sgr dSph without dark matter are presented in section|3] Some of 
these galaxies were simulated in different contents of dark matter. These results are shown in Section B] Finally we 
summarize some features of both the progenitor of Sgr dSph and its initial dark halo. 

2. Model 

2.1. Initial conditions 

There have been several efforts to unravel the real shape of the dark matter halo of The Milky Way. Some studies 
done by Ibata et al. Il23l and Majewski et al. [36] suggest that its shape must be close to spherical; Helmi ifTTl argues 
that the velocities of the leading tidal tail favor a prolate dark matter halo; nevertheless, the over-densities located at 
90 kpc from the galactic center found by Newberg et al. Il44l suggest that the dark halo is nearly spherical or slightly 
oblate ifBTl . The mystery seems to be solved taking into account the bifurcations in the tidal tails observed by Be- 
lokurov et al. [4| and Yanny et al. BTI . Fellhauer et al. [ 13 1 found that, if the potential of the dark halo were exactly 
spherical, the debris of Sgr dSph would lie in a single plane and no bifurcation would exist; furthermore, If the halo 
is too oblate or too prolate, then debris are scattered over a wide range of locations. Furthermore, if the simulated 
satellite has a considerable mass (more or equal to 7.5 x 10 8 M ), the bifurcations blurs. Penarrubia et al. PTll found 
that bifurcations in the leading tail of Sgr dSph stream naturally appear in models where the dwarf is represented 
by an exponential disk and if it is misaligned with respect to the orbital plane; furthermore, Law and Majewski [35] 
reproduced the majority of major constraints of the angular positions, distances and radial velocities of the young tidal 
debris streams, introducing a non axisymmetric potential in the form of a nearly oblate dark matter halo. With these 
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Figure 1 : Rotation curve of the Milky Way using the parameters of Johnston et al. 1271 . 



results, both studies make a realistic picture of Sgr dSph. 

Since we aim to find out if Sgr dSph must have dark matter for its evolution and what was its initial content, we 
will only focus on reproducing the physical properties observed in the main body, not in reproducing the bifurcations 
of the tidal tails. To do that, the Milky Way is modeled as a three-component rigid potential, in which the disk is 
represented by a Miyamoto & Nagai potential Il42ll . the spheroid by a Hemquist potential 1 18 1 and the dark halo by a 
spherical logarithmic potential 1 28 1 : 
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Where the disk mass M c i = 1 x 10" M Q ; the spheroid mass M s = 3.4 x 10 1() M ; the dark halo circular velocity 
v/, = 128 km/s; and the lenght scales related to each components of the Milky Way: a=6.5 kpc; b=0.26 kpc; c=0.7 
kpc; d=12 kpc. These parameters provide a nearly fiat rotation curve between 1 and 30 kpc and a disk scale height of 
0.2 kpc for our galaxy as can be seen in figure [T] 

If the Milky Way is modeled as a rigid potential, we can isolate the effects of tidal interactions on the internal 
structure and kinematics of the dwarf 113T1 ; although, with this approach, we neglect the effects of dynamical friction, 
the evolution of the potential of the host galaxy and interactions between dwarfs. 

With the purpose of establishing the initial conditions that could have the progenitor of Sgr dSph, a test particle 
was evolved under the action of the rigid potential, with the current position and radial velocity of Sgr dSph; that is 
16 kpc and 171 km/s respectively. Since the tangential velocity can not be measured, it is taken as a free parameter; 
we varied it from -325 to 325 km/s because with these values the apocenter of Sgr dSph will be between 40 and 60 
kpc, which is in agreement with the tidal tails found by Martmez-Delgado et al. Il38l . Using numerical simulations, 
Martmez-Delgado et al. [38] found that these tidal tails were disrupted from the satellite approximately 6 Gyrs ago. 
This period of time is very advantageous since some galactic evolution and dynamical friction effects can be neglected. 
EH. 

Higher values in the tangential velocity will cause the satellite to move farther from the galactic center; hence 
dynamical friction effects for the evolution of Sgr dSph should be taken into account. Thus, we will assume that Sgr 
dSph has been orbiting the Milky Way for a period of time smaller than 7 Gyr. 
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Table 2: Initial conditions that reproduce, after a time, the 
actual galactocentric position and radial velocity of Sgr 
dSph. A is the apocenter; P the pericenter; V is the tan- 
gential velocity of the particle at the apocenter and eg, the 
excentricity. 

A (kpc) P (kpc) V(km/s) io~~ 



44 12 


-104 


0.56 


45 12 


-102 


0.56 


46 12 


-101 


0.56 


47 12 


-100 


0.57 


48 12 


97 


0.59 


49 12 


96 


0.59 


50 12 


95 


0.60 


51 12 


94 


0.60 


52 12 


93 


0.61 


53 13 


-90 


0.60 


54 13 


-90 


0.60 


60 13 


-86 


0.63 



For each tangential velocity value, we obtained a Rosette-like orbit for the motion of the test particle; from each 
path, we computed its apogalactic distance, the tangential velocity at this point and its pericenter distance. Each of 
these results could become a set of possible initial conditions for the progenitor of Sgr dSph. In order to establish 
which of these values reproduce its current galactocentric position and radial velocity, the target particle was located 
initially at each apogalactic distance in the Y axis with an initial velocity given by the velocity at that point. This 
particle was evolved under the influence of the rigid potential. The values which reproduce the current position and 
radial velocity of Sgr dSph, are shown in table [2]. 

As we can see in table [5] these results are in good agreement with studies made by Ing-Guey and Binney 1 26 Q 
Ibata and Lewis [21 1; Law et al. ||32ll33l and Majewski et al. 1571 . Sgr dSph apocenter is between 44 and 60 kpc and 
its pericenter, between 12 and 13 kpc. The orbit period is 0.7 Gyrs. We also found that the present tangential velocity 
of the dwarf should be 298 + 27 km/s. 



2.2. Dynamical friction 

Since we modeled the Milky Way as a three component rigid potential, we can not simulate massive systems. We 
used the equation of Binney and Tremaine [5]|^] To estimate the maximum mass that we can use in our simulations: 

/2.64xlO I1 \/ r \ 2 1 v c \/10V\ 

M =h^HM [^^)[i^rr (4) 

Where Coulomb logarithm is log A=10 for the case when a globular cluster is decaying into the center of the 
Milky Way, r and v c are the initial position and the circular velocity at this point of the satellite respectiveljj^] For our 
simulations we take r=44-60 kpc. From our model of the rigid potential, the circular speed for each position is v c = 
207.802 - 202.345 km/s respectively. 



In this study they found two consistent solutions for the Srg past: one which it has an initial dark halo that caused to sink inwards under the 
influence of dynamical friction and another where Sgr was moved at all times on the same short period orbit. 

-Despite of this equation is valid only when the host galaxy is modeled as an isothermal sphere, we will use it because most of the orbit of Sgr 
dSph lies in the galactic halo. 

3 This calculation assumes the cluster maintains a circular velocity as it spirals to the center of the Milky Way. 
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Since our simulation time is f„ m * 7 Gyr, the falling time or tf a u should be much higher, we chose T fail = lOGyr. 

With those data, we found two maximum values for the initial mass of the system. Ifr = 44kpc, M max = 1.06xl0 9 
M ; when r = 60 kpc, M max = 1.92 X 10 9 M Q . 

2.3. Satellites and its dark halo 

Although the representation of Sgr dSph as a disk, reproduces a realistic picture for its evolution, there is no 
evidence yet of significant rotation in the main body of the satellite [35 1; that is why the progenitor of Sgr dSph was 
modeled as a Plummer sphere |48 1 whose potential and density are described respectively by: 
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where M sat is the initial satellite mass and ro its Plummer radius. We used the Aarseth algorithm [1,2] in order to 
construct the Plummer sphere. The initial mass of the satellites was varied from 5 x 10 8 to (1 - 1.9) x 10 9 M , since 
the present mass of Sgr dSph is (2.8 - 3.7)xl0 8 M Q Il33ll . The upper limit for the initial mass has been chosen to avoid 
dynamical friction effects as it was shown in the last section. The Plummer radii used were 0.3,0.5,0.6 and 1.2 kpc. 

In the case of satellites with a dark matter component, the dark matter halo was modeled using the method of 
distribution functions proposed by Mashchenko and Sills [40 1 which creates a dark halo as a spherical distribution 
of particles that describes a Navarro-Frenk & White (NFW) density distribution. The stellar component for these 
systems was modelled as a Plummer sphere. Since this method provides a distribution function consistent only with 
dark matter, we virialized the system composed by the Plummer sphere and the NFW halo with Gadget2 [49 1 until it 
reaches stability. Again, the maximum initial mass used for simulating systems with barionic and dark matter compo- 
nents was (1 - 1.9) x 10 9 M Q to avoid dynamical friction effects. 

The simulations of the evolution of the satellites with only barionic matter and systems of stellar and dark matter 
components, were performed using a modification of the Gadget-2 code [49]. On it, we inserted the rigid potential 



described in section 2. 1 to model our galaxy. All the systems were located initially at 44 and 60 kpc from the galactic 
center, with velocities as given in table [2] The softenings lenght used in our simulations were: 0.015 and 0.01 for the 
halo and stars respectively. 



3. Galaxies without dark matter 

All the satellites with only barionic matter were virialized during 1 Gyr using the Gadget-2 code in order to stabi- 
lize the systems; afterwards, we located them at 44 and 60 kpc from the galactic center with the respective velocities 
given in table [2] Despite of we are not taking into account dynamical friction effects, the current position and radial 
velocity of Sgr dSph are reproduced in our N body simulations. 

The current physical properties of Sgr dSph shown in table [T] were compared with the same quantities calculated 
in our satellites in order to find out the features of the progenitor of this dwarf galaxy. Theoretical half light radius, 
velocity dispersion, surface brightness and mass to light ratio are determined by the method shown by Klessen and 
Kroupa 1301 . in which these parameters are calculated by a terrestrial observer. Tables [3] and |4] show the results of 
our simulations when the initial position of the dwarfs is 44 and 60 kpc respectively. For both tables, t corresponds 
to the time at which the simulated satellites are located at the current position of Sgr dSph and have its observational 
galactocentric radial velocity; Mass is the mass of these systems at time t. ri/2, fiQ, cr and M/L are the half brightness 
radius, surface brightness, velocity dispersion and mass to light ratio respectively measured at time t as well. 
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Table 3: Physical properties of simulated satellites when they are located initially at 44 kpc from the galactic center. Both 
Initial mass, M sa , and mass obtained at time t, are given in M Q ; the half light radius ri/2 and Plummer radius rg are in kpc; 
time when Sgr dSph is located at the current galactocentric position f, is in Gyr; the central surface brightness /jo is measured 
in L e /kpc 2 ; the units of velocity dispersion and mass to light ratio are km/s and M G /L Q respectively. 



M sat 


ro 


t 


Mass 


n/2 






M/L 


1 x 10 y 


0.6 


0.9 


1 x 10 9 


0.18 


8.18 x 10 8 


33.17 


2.45 


9 x 10 8 


0.6 


0.9 


7.2 x 10 8 


0.18 


7.24 x 10 8 


31.16 


2.43 


8 x 10 8 


0.6 


0.9 


6.4 x 10 8 


0.18 


6.51 x 10 8 


29.08 


2.38 


7 x 10 8 


0.5 


0.9 


7 x 10 8 


0.15 


7.99 x 10 8 


30.39 


2.43 


6x 10 8 


0.5 


0.9 


6x 10 8 


0.15 


6.85 x 10 8 


28.16 


2.44 


5 x 10 8 


0.5 


0.9 


4x 10 8 


0.15 


5.81 x 10 8 


25.39 


2.38 


5x 10 8 


0.3 


0.9 


5 x 10 8 


0.10 


51.54 x 10 9 


34.19 


2.39 


7x 10 8 


0.6 


0.9 


5.6 x 10 8 


0.18 


5.33 x 10 8 


26.92 


2.43 


6x 10 8 


0.6 


0.9 


4.8 x 10 6 


0.18 


4.65 x 10 8 


24.74 


2.39 


5 x 10 8 


0.6 


0.9 


5 x 10 8 


0.17 


4.07 x 10 8 


22.17 


2.29 


1 x 10 9 


1.2 


0.9 


4x 10 8 


0.32 


1.30 x 10 8 


21.65 


3.71 


7 x 10 8 


1.2 


0.9 


2.1 x 10 8 


0.34 


5.8 x 10 7 


17.02 


4.8 



Table 4: Physical properties of simulated satellites when they are located initially at 60 kpc from the galactic center. The 
notation is the same as table[3] 



M sat 


10 


t 


Mass 


ri/2 


HO 


0"o 


M/L 


1 x 10 9 


0.6 


1.2 


1 x 10 9 


0.18 


7.8 x 10 8 


33.69 


2.58 


6x 10 8 


0.6 


1.2 


4.8 x 10 8 


0.19 


4.3 x 10 8 


25.11 


2.53 


6x 10 8 


0.5 


1.2 


6x 10 8 


0.16 


6.63 x 10 8 


28.5 


2.62 






3.5 


4.8 x 10 8 


0.15 


6.61 x 10 8 


27.05 


2.41 


5 x 10 8 


0.6 


1.2 


4x 10 8 


0.18 


3.81 x 10 8 


22.44 


2.38 


5 x 10 8 


0.5 


1.2 


5 x 10 8 


0.15 


5.56 x 10 8 


25.6 


2.47 






3.5 


3.5 x 10 8 


0.14 


5.49 x 10 8 


23.96 


2.32 
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Figure 2: Lagrange radii of some simulated satellites when they are located initially at 44 kpc from the galactic center. In all graphics, lines from 
bottom to top, represent the radius that contains 1%, 10%, until 90% of the initial mass. When the radius which represents 1% of the initial mass 
of the galaxy reach the top of each graphic, the satellite is totally disrupted. The maximum number shown in the Y axis, represents the initial virial 
radius of the satellites, which is calculated using the expression R = r u C where r u is the Plummer radius and C=5 is the cutoff taken to construct 
the Plummer spheres. 

Comparing the current observed properties of Sgr dSph with the results from our sample of simulated galaxies, 
we can see in tables [3] and [4] that most of them have a higher mass at time f; except from the extended satellites (those 
which have r ( > = 1.2 kpc) and the ones with initial mass of 5 x 10 8 M . This is not surprising since the satellites 
reproduce the current galactocentric position and radial velocity of Sgr dSph at early stages of their evolution. 

The half light radius, r\j%, is lower for all the satellites except again for the extended ones (ro = 1.2 kpc). These 
satellites reproduce the measured r\j% of Sgr dSph at 3 Gyrs. Nevertheless, at this time they do not reproduce the 
galactocentric position and radial velocity of Sgr dSph. 

The surface brightness and the velocity dispersion are larger than the ones measured for Sgr dSph. 

Figure|2]shows the lifetimes for some of the simulated satellites when they are located at 44 kpc from the galactic 
center. For the case of satellites whose initial mass is very high (1 x 10 9 ) M or its Plummer radius is very small 
(0.3 kpc) more than 30% of their initial mass survive for 10 Gyrs. It means that they will live for a Hubble time; 
furthermore, from the same figure, we can see that these galaxies, between 6 and 7 Gyrs, have masses higher than 
the current accepted value for Sgr dSph. It means that these satellites are more stable than the real Sgr dSph. This 
behavior is present in the first galaxies listed in table [3] with this evidence, we conclude that the progenitor of Sgr 
dSph could not have been a very dense galaxy. In the case of a satellite of initial mass of 6 x 10 8 M Q and Plummer 
radius of 0.6 kpc, this galaxy is totally disrupted before 10 Gyrs. The same behavior is presented in satellites with 
initial mass of 5 x 10 8 M Q and Plummer radius of 0.6 kpc. At 6 Gyr both galaxies reproduce the current mass of Sgr 
dSph; nevertheless, they fail in reproducing the other properties of Sgr dSph at this time. 

Figure [2] also shows that satellites with Plummer radius of 1.2 kpc are destroyed before 5 Gyrs, as we can see in 
the same figure; but these systems are the only ones which reproduce the half brightness radius of Sgr dSph. 

If the satellites are initially located at 60 kpc from the galactic center, their lifetimes are much larger, as we can 
see in figure [3] where we show the evolution of the Lagrange radii that contain from 1% to 90% of the initial mass of 
a galaxy with initial mass of 6 x 10 8 M Q and a Plummer radius of 0.6 kpc. This galaxy is totally disrupted before 10 
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Figure 3: Evolution of a satellite with initial mass and Plummer radius of 6 X 10 M Q and 0.6 kpc respectively. This satellite initially is located at 
a galactocentric distance of: Top: 44 kpc. Bottom: 60 kpc. 

Gyrs if it is located at 44 kpc; nevertheless, it can survive for a Hubble time if its initial position is 60 kpc, because 
30% of its initial mass survives for 10 Gyrs. The same behavior is observed in most of the galaxies in our sample, 
except for satellites with Plummer radius of 1.2 kpc. Their lifetimes are always lower than 5 Gyrs whenever their 
initial positions are. 

Klessen and Kroupa 1 30 1 ; Casas and Kroupa |6| and Casas et al. |7| have argued against the existence of dark 
matter in dwarf galaxies. In those studies, they obtained large values in the M/L ratio for spheroidal systems without 
a dark component. This can be explained if it is assumed that dwarf galaxies are out of equilibrium systems. In this 
study, we found a remarkable different result: the M/L ratio for all the simulated satellites without dark matter, is 
lower than the value measured for Sgr dSph. This could be explained if this satellite has a different origin than the 
other dwarf galaxies of the Milky Way. 

On the other hand, many observational studies have failed trying to find out evidence of dark matter in dwarf 
galaxies |[T2l l9l. Nevertheless, none of our simulations without dark matter could reproduce the measured physical 
properties of Sgr dSph. This result suggest that this satellite, at the beginning of its evolution, should had been 
immersed in a dark matter halo. Furthermore, the simulations showed in this section point out that its barionic 
component could not have been a very dense object. 



4. Galaxies immersed in a dark matter halo 



As we mentioned in section 2.3 systems composed by barionic and dark matter, were virialized under their self 



gravity using Gadget 2 [49] until they became stable. Figure |4] shows this virialization during 8 Gyrs. 

The concentration parameter in dark matter halos depends on its mass Il34ll40l . redshift [40 1 and the initial power 
spectrum of density fluctuations 1341 . The mass concentration (M-c) relation of virialized halos has been determined 
by various authors using N-body simulations IfTD or observational data 



Using the M-c relation derived by Duffy et al. ifTTl . Comerford and Natarajan [8] and Oguri et al. 114611 . we com- 
puted the concentration parameter of the dark matter halo that could have had the progenitor of Sgr dSph 7 Gyrs ago. 
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Figure 4: Relaxation of a system with a dark halo of 8 X 10 M Q , r s = 3 kpc and concentration of 10. The barionic component is a satellite of 
7 X 10 s M e and Plummet radius of 0.6 kpc. Note the collapse of the system during the first 3 Gyr. The system reaches stability after 7 Gyrs. 



We used 5 x 10 8 M and 1.9 x 10 9 M Q for the calculation^] Using the M-c relation of Duffy et al. ifTTl . we obtained 
c= 4.52-4.06 for the value of masses respectively. Using the M-c relations of Comerford and Natarajan [8 1 and Oguri 
et al. [46 1 we obtained c= 12.36-10.12 and c= 11.72-10.51 respectively. Taking into account these values, we used 
for our simulations, concentrations parameters of 4 and 10; furthermore, since numerical simulations suggest that 
4 < c < 22 [34], we made one simulation using c=6. 

For the simulations of satellites with dark matter, we choose for the barionic component, two types of satellites: 
not very dense ( M < 7 x 10 s M Q ) and extended satellites (ro = 1.2 kpc). 

When we studied the dark component of our simulated systems, we observed that the disruption of a dark halo 
depends on both its initial mass and concentration parameter. As we can see in figure [5] in the case of a system whose 
dark component has a mass of 8 x 10 8 M and concentration of 10, the barionic component begins to be disrupted 
when the dark halo has lost 44% of its initial mass. At 2.2 Gyrs of evolution, only remains 10% of the initial mass 
of the dark matter halo; while the satellite solely has lost 20%. In the same figure we show the evolution of the same 
satellite only with barionic matter. As we expected, the dark component protects the galaxy against tidal disruption. 
When a system with a dark halo of initial mass of 8 x 10 8 M Q and concentration parameter of 4 loses 10% of its mass, 
the galaxy begins to be disrupted. At 2.8 Gyr of evolution, the dark halo contains only 1% of its initial mass while the 
stellar component contains 60%. 

We evolve each system immersed in a rigid potential for a time lower than 7 Gyrs; except for the systems with 
dark matter halos with 9 x 10 6 ~ 7 M Q and barionic component of 7 X 10 8 M Q , because the fall time given by equation 
is larger than 10 Gyr when they are located initially at 44 kpc from the galactic center. The half light radius, surface 
brightness, velocity dispersion and mass to light ratio of two galaxies with initial mass of 7 x 10 8 M and Plummer 
radii of 0.6 and 1.2 kpc are presented in table|5] As in table|5] those parameters correspond to the physical properties 
observed by a terrestrial observer when each system reproduces the current galactocentric position and radial velocity 
of Sgr dSph. Continuous lines represent total disruption of the simulated satellites. 

Simulated galaxies with Plummer radii of 0.6 kpc have a mass consistent with the current range measured in Sgr 
dSph, as we can see in table|5]after 5.7 and 6.5 Gyr of evolution. Nevertheless, their mass to light ratio is lower than 
the measured value in Sgr dSph ; while their surface brightness is higher (these satellites are more luminous than Sgr 
dSph). Only when galaxies are immersed in less massive dark halos, i.e. 9 x 10 7 or 9 x 10 6 M Q , they reproduce most 
of the physical properties of Sgr dSph. This result suggests that its initial content of dark matter could have this range 



4 Those are the minimum and maximum values used in the simulations. 
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Table 5: Physical properties of a galaxy with an initial mass of 7 X 10 8 M Q and Plummer radius of 0.6 and 1.2 kpc when is 
immersed in different contents of dark matter. Its initial position is 44 kpc from the galactic centre. Columns have the same 
meaning as table 3; furthermore, the units of the physical parameters are equal. Note that the half brightness radius of the 
satellites does not vary with the inclusion of a dark matter halo. 



ro 




Halo 




t 


Mass 




Mo 


O"0 


M/L 




M[M Q ] 


r s [kpc] 


C 
















8 x 10 8 


3 


10 


5.7 


2.1 x 10 8 


0.14 


4.71 x 10 8 


20.20 


1.95 




8 x 10 8 


7.5 


4 


0.9 


5.6 x 10 8 


0.18 


6.36 x 10 8 


28.77 


2.30 










6.5 


2.8 x 10 8 


0.16 


4.1 x 10 8 


20.31 


2.05 


0.6 


9x 10 7 


7.5 


4 


0.9 


5.6 x 10 8 


0.19 


6.53 x 10 7 


27.31 


16.56 










8.2 


1.4 x 10 8 


0.16 


2.72 x 10 7 


15 


16.9 




9x 10 6 


7.5 


4 


0.9 


5.6 x 10 8 


0.18 


7.18 x 10 6 


27.27 


184 










7.7 


2.1 x 10 8 


0.16 


3.68 x 10 6 


16.9 


160 




1 x 10 7 


7.5 


4 


2.6 
5.7 


7x 10 6 


0.65 


1.9 x 10 4 


21.84 


22066 


1.2 


8 x 10 7 


8.3 


6 


4.4 


7x 10 6 


0.69 


2.2 x 10 3 


167 


600 




8 x 10 8 


6 


10 


4.4 


7 x 10 6 


0.39 


7.8 x 10 s 


93.87 


500 
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Table 6: Physical properties of a galaxy with an initial mass of 1 X 10' M Q and Plummer radius of 1.2 kpc when is immersed 
in different contents of dark matter. Its initial position is 44 kpc from the galactic centre. Third column shows the mass of the 
barionic matter at time t. 





Halo 




t 


Mass 


1/2 


Mo 


O"0 


M/L 


M[M ] 
9x 10 7 


r s [kpc] 
7.5 


C 

4 


0.9 
5.9 


1 x 10 8 


0.45 


5.8 x 10 5 


17.86 


406 


5x 10" 


6 


10 


0.9 


5 x 10 s 


0.3 


9.4 x 10' 


22.03 


6 



of masses. Nevertheless, when we evolve satellites with initial dark matter halos of 1 x 10 7 and 8 x 10 7 M and the 
barionic component with a Plummer radius higher than 0.6 kpc, we can see in table[5]that they can not survive for a 
time larger than 5 Gyrs; that is, these simulated galaxies are less stable than Sgr dSph. This result rejects the fact that 
the initial content of dark matter of Sgr dSph could have been of the order of 10 7 or 10 6 M Q . 

All of the simulated satellites with Plummer radius of 0.6 kpc survive for 10 Gyrs and as we can see in Figure [5] 
if the satellite is immersed in a dark matter halo of concentration 10 and initial mass of 8 x 10 8 M Q , at least 10% of 
their initial mass will survive for a Hubble time. 

When we compare the half brightness radius of satellites with Plummer radius of 0.6 kpc (table |4]> and 1.2 kpc 
(tables|4]and[5]) we can see that only these galaxies reach a value similar to the observed in Sgr dSph in some of their 
stages of evolution. This suggests that at the beginning of the evolution of this dwarf galaxy, it should have been an 
extended satellite with a Plummer radius larger than 0.5 or 0.6 kpc. 

The evolution of extended satellites (Plummer radius of 1.2 kpc) immersed in a dark matter halo is presented in 
figure [6] As we noted before, their lifetimes are not higher than 5 Gyrs even when their halos are massive and have 
a concentration parameter of 10. Since the radii of these dark matter halos can not be larger, because they could 
produce some notable effects in the Milky Way, both their initial masses and concentration parameter must be larger 
than the values used here. It implies that in order to figure out what was the initial dark matter content of Sgr dSph 
using N Body simulations, dynamical friction effects should be taken into account. The mass of those halos, could 
possibly had an order of 10 10 or 10 11 M Q . In one of his recent papers, Nierdeste-Ostholt et al. [45 1 conclude through 
measurements of the total luminosity of the main body and tidal tails of Sgr dSph that the initial mass of its dark 
matter halo was su 6.3 x 10 9 M Q . This result could be tested using N Body simulations to see if one extended galaxy 
with this halo can reproduce all the measured physical properties in Sgr dSph. 



5. Conclusions 

We have shown several N body simulations of the evolution of the Sgr dSph with and without dark matter. In the 
simulations of satellites with only barionic matter, we found that very dense satellites are more stable than Sgr dSph 
because during the entire simulation time, that is, 7 Gyrs, their masses are always larger than the value measured in 
Sgr dSph. On the other hand, galaxies with initial masses of (5 - 6) x 10 8 M Q reproduced the current bound mass 
of Sgr dSph at 6 Gyrs. Extended satellites (ro ~ 1.2 kpc) were the only able to reproduce the half light radius of 
Sgr dSph. Despite these results, it was not possible to obtain the entire properties observed in Sgr dSph by evolving 
systems with only a barionic component; that is why we conclude that the progenitor of Sgr dSph should had been 
immersed in a dark matter halo, as well as the current satellite; despite the fact that many observational studies have 
failed in finding observational evidence of dark matter in dwarf galaxies. 

Simulations of systems with dark and barionic matter suggest that the progenitor of Sgr dSph had to be an extended 
galaxy (ro > 0.6 kpc) in order to reproduce the observed half light mass radius; furthermore, the initial mass of its 
dark matter halo had to be larger than 10 8 M Q . 
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2 

Time [Gyr] 



4 



2 4 

Time [Gyr] 



Figure 6: Evolution of satellites with Plummer radius of 1.2 kpc and initial mass of: left: 7 X 10 s M e . Right: 1 X 10 9 M Q . Note that there are no 
difference between a satellite with only barionic matter and the same one immersed in a dark halo. The satellite of 7 X 10 8 M Q has a dark halo of 
8 X 10 8 M while the satellite of 1 X 10 9 M Q , another one with 5 X 10 8 M Q . Both have a concentration parameter of 10. Note that these are systems 
with the maximum allowed mass: 1.5 X 10 9 M Q . 
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